Chemical versus van der Waals Interaction: The Role of the Heteroatom 
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We performed first-principles calculations aimed to investigate the role of an heteroatom like 
N in the chemical and the long-range van der Waals (vdW) interactions for a flat adsorption of 
several tt— conjugated molecules on the Cu(110) surface. Our study reveals that the alignment of 
the molecular orbit als at adsorbate-substrate interface depends on the number of heteroatoms. As 
a direct consequence, the molecule-surface vdW interactions involve not only 7r-like orbitals which 
are perpendicular to the molecular plane but also cr-like orbitals delocalized in the molecular plane. 

PACS numbers: 73.20.-r, 68.43.Bc,71.15.Mb 



In the last decade, the molecular electronics emerged as 
a promising field able to provide the technology necessary 
to develop devices such as organic light-emitting diodes 
(OLEDs) 0, Q, organic field-effect transistors [3, 0, [H( 
or ultra-high-density memory circuits p. In this con- 
text, the ability to reliably describe the electronic prop- 
erties 0, | or the mechanical manipulation process of 
molecules 0, [To[ on surfaces is essential to understand 
and design the functionality of such devices. The achieve- 
ment of these goals strongly depends on the accuracy of 
the state-of-the-art theoretical methods used to assess 
the interaction between a molecule or a molecular layer 
and a substrate of choice. Nowadays, the density func- 
tional theory (DFT) is the theoretical tool of choice to 
analyze [7|, |S|, LUJ, LL2[ and to predict [13, LL4| the electronic 
properties of systems characterized by strong chemical 
bonds. 

However, the loosely bounded physical systems 



15] 



represent one of the major challenges for DFT because 
the effective Kohn-Sham (KS) potential does not exhibit 
the correct asymptotic behavior. In particular, the cur- 
rently used exchange correlation energy functionals like 
the local density approximation (LDA) or the generalized 
gradient approximation (GGA) do not properly describe 
the long-range van der Waals (vdW) interactions. For in- 
stance, the GGA fails to predict a bonding ground-state 
for a van der Waals system like graphite [161 ]. 

One way to circumvent this limit is to include the 
dispersion effects in DFT in a semiempirical fashion, in 
which the total energy of the physical system is a sum of 
the self-consistent Kohn-Sham energy as obtained from 
the DFT calculations and semiempirical dispersion cor- 
rection which depends on the interatomic distance. For 
example, a semiempirical treatment of the vdW disper- 
sions in ab initio calculations was used to investigate, for 
instance, the adenine on graphite(OOOl) [17]. 

Another way to accurately account for the vdW inter- 
actions is a first-principle approach [l8[ which requires 
the calculation of a nonlocal correlation energy functional 



as a post GGA perturbation procedure. This method 
was successfully employed to investigate the benzene, 
naphthalene and phenol on graphite (0001) [l9[ and a- 
Al 2 O 3 (0001) or thiophene on Cu(110) surfaces. 

In this Letter, we investigate the role of the het- 
eroatoms on the chemical and the van der Waals interac- 
tions for a fiat adsorption geometry on the Cu(110) sur- 
face of three prototype 7r— electron systems as benzene 
(Bz), pyridine (Py) and pyrazine (Pz) molecules. Our ab 
initio simulations reveal that, for the Bz the long-range 
dispersion effects are basically important only for the ad- 
sorption energy. This is not the case of the Py and Pz, 
where due to their low lying tt— orbitals, the inclusion of 
the long-range correlations drastically influence the ad- 
sorption geometry and electronic structure, i. e. the Py 
becomes chemisorbed on the surface while the Pz binds 
to the surface mostly due to the van der Waals interac- 
tion. To evaluate the dispersion effects we used the ab 
initio (vdW-DF) [l8[ as well as a semiempirical (DFT-D) 
21, method. In particular, a key result of our study is 



that a state of the art investigation of the bonding mecha- 
nism leading to a flat molecule— substrate adsorption ge- 
ometry requires the use of both methods to correctly de- 
scribe the adsorption-geometry, the electronic structure 
and the correlation effects of the molecule— surface inter- 
face. Moreover, we show that the nature of the molecular 
orbitals involved in the dispersion interactions depends 
on the alignment of these orbitals with respect to the 
Fermi level of the adsorbate— surface system. 

Our first— principles total-energy calculations are car- 
ried out in the framework of density functional theory 
(DFT) by employing the generalized gradient approxi- 
mation (GGA) in a pseudopotential plane- wave formu- 
lation (projector augmented wave method) [23[ as im- 
plemented in the VASP code \2% l25j |. In particular, in 
our study we employed the PBE [26[ flavor of GGA. The 
plane— wave basis set consists of all plane waves up to a 
kinetic energy of 450 eV. The molecule— Cu(110) system 
is modeled within the supercell approach and contains 
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FIG. 1: (Color online) Energetic levels diagram for Bz (left), 
Py (middle) and Pz (right) molecules. Due to the presence 
of N atoms in the Py and Pz molecules the tt— orbitals are 
shifted to lower energies while some of a— orbitals are pushed 
to higher energies [341 ]. 



five atomic layers of copper with the adsorbed molecule 
on one side of the slab [27[. It was generated with the 
theoretical bulk copper lattice parameter of 3.63 A and 
has a p(4x6) in-plane surface unit cell. During our ab 
initio calculations, the uppermost two copper layers and 
the molecule atoms were allowed to relax until the atomic 
forces are lower than 0.001 eV/A. 

The evaluation of the vdW forces by using the DFT- 
D 21, [22[ method in the self-consistent cycle of our ab 



initio calculations are mandatory to obtain the proper 
equilibrium adsorption geometry and the corresponding 
electronic structure of the Py- and Pz-Cu(llO) systems. 
Moreover, we also evaluated the correlation effects by 
using the vdW-DF functional [l8[ as a post processing 
procedure (i.e., non-self-consistent) on a charge density 
obtained from the standard GGA calculations [28] using 
the [IK code developed in our group. 

As compared with Bz (CeHe), in the Py (C5H5N) 
molecule a CH group of atoms is replaced by a TV atom, 
while the Pz (C4H4N2) has two CH groups replaced by 
two N atoms. All molecules contain a 6tt— electron sys- 
tem and they can adsorb with the molecular plane par- 
allel to the Cu(110) surface 30l 3lL 32]. From electronic 
point of view, it is well known [33|, l34| that the presence 
of TV atoms in the conjugated heterocycle molecules lower 
their tt— orbitals while some of a orbitals are pushed to 
higher energies. This peculiar feature is clearly shown 
in the energy levels diagram of the isolated Bz, Py, Pz 
molecules (see Fig. [1]). 

The Bz is considered to be chemisorbed on the Cu(110) 
surface [13, HH due to an effective hybridization of the 



p z — atomic orbitals (those perpendicular to the molec- 
ular plane) with the Cu d— bands. Note that these 
p z — atomic orbitals are forming the tt— orbitals of the iso- 
lated molecule. On the contrary, in the case of Py and Pz 
molecules the presence of N atoms lower their tt— orbitals 
which implies that the interaction of the corresponding 
p z — atomic orbitals with the Cu d— bands is expected to 
be much weaker [33[. Indeed, using only the standard 
GGA approximations (PBE) in our DFT calculations of 
the flat adsorption geometry of Py and Pz molecules on 
Cu(110), we find that the adsorption site is not sensitive 
to the underlying geometry of the surface. The average 
Py(Pz)— surface metal distance is 2.98(2.99) A, while for 
the Bz this distance is 2.43 A. The magnitude of Py- 
and Pz-surface separation distances suggest that these 
molecules are physisorbed on Cu(110), while as already 
mentioned, the Bz is chemisorbed. 




FIG. 2: (Color online) Side views of the relaxed geometries 
when using standard DFT (left panel) and DFT+D (right 
panel). The inclusion of the van der Waals interactions at a 
semi-empirical level in the DFT+D approach [22|] decreases 
significantly the molecule— surface separation distance in the 
case of the Py and Pz molecules. The adsorbate-substrate 
distances are given in A. 

In the next step we performed geometrical relaxations 
of the molecule— metal systems including the van der 
Waals forces using the DFT-D approach. Since the vdW 
forces are attractive, their effect is to bring the molecule 
closer to the surface (see Fig. [2]). For the Bz molecule 
this effect on the adsorption geometry is not significant, 
since the average Bz— surface distance is decreased from 
2.43 to 2.35 A. On the contrary, the inclusion of the long 
range dispersion effects has a huge impact on the geom- 
etry of the Py and Pz molecules adsorbed on Cu(110) 
surface. The average Py(Pz)-surface distance decreases 
significantly from 2.98(2.99) to 2.43(2.59) A. Also the tilt 
angle of the Py molecular plane changes from +3° (DFT) 
to w -7° (DFT+D). 

This structural change drastically affects the electronic 
structure of the Py- and Pz-Cu(llO) systems. In Fig. [3] 
we present the local density of states (LDOS) for Bz, Py 
and Pz molecules in the case of the relaxed geometries 
before (DFT) and after inclusion of the van der Waals 
forces (DFT+vdW). At the molecular site the basic char- 
acteristic displayed by the Bz— Cu(110) system are the 
broad bands with ^1^,3— character, while the a ones are 
quite sharp and very localized on the molecule. Also 
this electronic structure is not affected by the inclusion 
of the long-range dispersion effects because the geometry 
of the system changes only slightly. On the contrary, for 
the Py and Pz molecules, due to a large molecule— metal 
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FIG. 3: (Color online) Upper panel: Local density of states (LDOS) of the Bz, Py and Pz molecules adsorbed on Cu(llO). In 
the case of Bz molecule, the electronic structure of the molecule— surface system is not significantly modified by the inclusion of 
the dispersion effects. This is not the case of the Py- and Pz-Cu(llO) systems, where the van der Waals attractive interaction 
brings the molecule closer to the surface and allows a stronger hybridization between the p z — orbitals and the Cu d— bands. 
Lower panel: Binding energies due to nonlocal correlation effects evaluated with the vdW-DF method. Red color indicates a 
strong non-local (NL) contribution to the binding energy while green denotes a much weaker contribution. In the case of the 
Pz molecule this contribution is delocalized over the whole heterocycle ring while in the case of the Bz and Py it mainly arises 
from the two carbon sites that directly bind to the surface. 



distance in the DFT relaxed geometry, the tt bands are 
less broad being localized mainly on the molecules (see 
Fig. [3]). The inclusion of the attractive van der Waals 
interaction brings the Py or Pz closer to the surface and 
thus it allows the p z — orbitals to hybridize more strongly 
with the Cu d— bands. This molecule— surface interaction 
results in broad bands with mixed 7r and metallic char- 
acter as already seen for the Bz— Cu(llO) system. Note 
also the appearance in the case of Pz of a broad g\ type 
band. 

To quantify the importance of the van der Waals in- 
teractions in the case of the flat adsorption process of 
the Bz, Py and Pz molecules on the Cu(llO) surface, 
in Table [J we report the calculated interaction E[ nt and 
adsorption E a( ^ s energies defined as: 

rp zp ( T^relaxed/ideal ^relaxed/clean \ 

Ant/ads — Ays ~ ^molecule ~+~ ^Cu(llO) r 

where E sys represents the total energy of the relaxed 
molecule— Cu(llO) system, E^£^ le represents the total 
energy of the isolated molecule and -^cu(iio) * s ^ ne ener gy 
of the Cu(llO) surface, both the molecule and surface be- 
ing in the same atomic configurations as in the relaxed 
molecule— Cu(l 10) system. Also, the ^molecule denotes 
the total energy of the equilibrium isolated molecule and 
-E'cufno) * s ^ ne ener &y °f the clean Cu(110) surface. 

The interaction energy EP£ t represents a measure of 
the strength of the chemical interaction between molecule 
and surface, while the adsorption one E^[ S T includes not 
only the energy gain due to the bond formation (EPf T ) 
but also the energy paid to deform the molecules and 



surface from their ideal configurations to those as found 
in the relaxed molecule— surface system. 





DFT relaxed 

7-iDFT t-iDFT 
Ant Ads 


DFT+vdW relaxed 

pDFT F DFT z-wciW-D pvdW-DF 
Ant Ads Ads Ads 


Bz 

Py 

Pz 


-301 -253 
-115 -108 
-35 -36 


-343 -252 -492 -543 
-196 -105 -417 -510 
+56 +70 -371 -577 



TABLE I: The interaction and adsorption energies of the Bz, 
Py and Pz molecules adsorbed on the Cu(110) surface (units: 
meV) for the geometries obtained using only DFT (left side) 
and DFT including the dispersion corrections (DFT+vdW, 
right side). Note that when the dispersion effects were con- 
sidered, the adsorption and interaction energies have been de- 
composed into the contribution arising from a DFT (E^ s T [nt ) 

and that given by the van der Waals corrections (^ads^nt D ^)' 

• pDFT+vdW _ pDFT , T ? vdW-D(F) 

Le -> Ads,int — Ads,int T Ads,int 

Although the molecular geometries are slightly dis- 
torted from a planar geometry after including the dis- 
persion effects as compared to the those obtained only 
with DFT, for Bz and Py the calculated adsorption en- 
ergies E^ S T are practically invariant Moreover, the in- 
teraction energies E^f T are lower for the DFT+vdW re- 
laxed geometries as compared to those calculated for the 
DFT relaxed ones. This theoretical finding implies that 
by including the van der Waals dispersions, the chemical 
interaction between the molecule and surface increasesl 
This effect is rather small for Bz since this molecule is al- 
ready chemically adsorbed on Cu(110) surface. However, 
the inclusion of the vdW interactions lowers significantly 
the interaction energy E^f T for Py on Cu(110). Also 
the average Py- and Bz-surface distances are now com- 
parable (~ 2.43 and 2.35 A, respectively, see also Fig. [2]). 



In consequence, the Py molecule becomes chemisorbed on 
the Cu(llO) surface. In this case, the the vdW attrac- 
tive forces are the key ingredient that simply triggers the 
chemisorption process of Py on this surface. As regard- 
ing the Pz molecule, the EP^ t indicates a weak bonding 
interaction (negative values) for the DFT relaxed geome- 
tries which turns into a repulsive interaction (positive 
values) after taking into account the dispersion effects. 
However, the adsorption energy due to the van der Waals 
interaction (i^ds^) * s mucn more negative as compared 
to EP^ T such that the total value of the adsorption en- 
ergy is negative. This means that the Pz molecule binds 
to the surface only through the van der Waals interaction. 

In order to deeply understand how the correlation ef- 
fects contribute to the molecule— metal interaction, in 
Fig. [3] we plotted the binding energies due to the non- 
local (NL) correlation effects as calculated within the 
vdW-DF theory [18[. For the Bz and Py molecules, the 
non-local contribution is mainly localized in the regions 
where the C atoms are situated on top of Cu atoms of the 
surface. In the case of Pz, the contribution of the non- 
local correlation effects is much more delocalized over the 
entire molecule. From the electronic structure point of 
view, this peculiar feature is related to a larger contri- 
bution to the NL part of the binding energy of the g\ 
molecular orbital in the case of the Pz with respect to 
the Py molecule. Note that for Bz this contribution is 
not significant since the o\ lies deep in energy with re- 
spect to the Fermi energy of the molecule-surface system. 
Therefore which molecular orbitals will contribute to the 
vdW interactions clearly depends on the specific align- 
ment of these orbitals at the molecule-surface interface. 
An animation is available as additional material [35|. 

To conclude, we prove that the van der Waals dis- 
persion effects together with the corresponding vdW 
attractive forces are crucial to reliably calculate in a 
self-consistent manner the proper equilibrium adsorp- 
tion geometry and the corresponding electronic struc- 
ture of 7T— conjugated heterocycle molecules adsorbed on 
Cu(llO) surface. The inclusion of the long-range dis- 
persion effects changes qualitatively the adsorption pro- 
cess of Py (C5H5N) on Cu(llO) surface from physisorp- 
tion to chemisorption, while for the Pz (C4H4N2) the 
van der Waals interactions represent the driving forces 
which bind the molecule to the surface. Our study also 
clearly pointed out that the vdW interactions depends 
on the alignment of the molecular orbitals at adsorbate- 
substrate interface and can involve not only 7r-like or- 
bitals (perpendicular to the molecular plane) but also 
a-like orbitals (in the molecular plane). Therefore we 
conclude that these effects will definitely play a key role 
in the adsorption process of other ir— conjugated hete- 
rocycle molecules as porphyrines or phthalocyanine on 
metal surfaces. 
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